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Abstract 


We exploit the relationship between multiobjective integer linear problem 
(MOILP) and data envelopment analysis (DEA) to develop an approach to 
a resource reallocation problem. The general purpose of the mathematical 
formulation of this multicriteria allocation model based on DEA is to enable 
decision-makers to take into account the efficiency of units under control 
to allocate additional resources for a new period of operation. We develop 
a formal approach based on DEA and MOILP to find the most preferred 
allocation plan taken account additional resources. The mathematical model 
is given, and we illustrate it with a numerical example. 
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1 Introduction 


The use of the data envelopment analysis (DEA) approach [1, 4] provides 
useful information for monitoring and management in a production orga- 
nization. The essential elements are knowledge reference units, sources of 
inefficient units, variations of productivity from one period to the next, and 
variable returns to scale. 

Multicriteria decision aiding (MCDA) and DEA are two techniques in 
operations research, which attract the attention of many researchers. Some 
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authors are interested in their relationship and how to exploit [2,3, 7-9, 11, 

13, 15,16]. We propose a mathematical model of multiobjective optimization 

based on DEA that can help managers or decision-makers adjust investments 

in different units of the control system by additional assignments to bring. 
The general purpose of the mathematical formulation of this multicriteria 

allocation model based on DEA is to enable decision-makers to take into ac- 

count the performance of units under control to allocate additional resources 

for a new period of operation. So, it is necessary to determine information 

intermediaries: 

— information on production or performance efficiencies; 

— information on the returns to scale presented by each unit; 

— optimal targets for each production unit. 


The mathematical model, we propose, takes into account all this infor- 
mation intermediaries. The objective is, on the basis of additional quan- 
tity of resources and demands expressed by the different DMUs, to allocate 
optimally and efficiently these resources available to maximize estimates of 
product returns based on information gathered from performance measures 
and returns of the ladder. In a given production system, the reference units 
(the efficient units) are those that use resources (inputs) rationally, at best, 
to achieve optimal results. This transformation of inputs into outputs is sup- 
posed to follow a function the best of which is that used by efficient units. 
These efficient units give an idea of best practices and thus provide the best 
possible production plans. On the basis of information on the units, namely 
the identification of the efficient and the inefficient, it is important to consider 
how this information can be used to assist the decision-maker in allocating 
resources for a new production period with the assumption that the system 
will use the same process of transformation. One of the information that 
can be provided by the application of the DEA method is the identification 
of the variable returns to scale associated with each unit of the system on 
a production function basis. Efficient and inefficient units can be either of 
increasing returns to scale, constant returns to scale, or decreasing returns 
to scale. 

In this paper, we present the modeling of the problem, the mathematical 
model and illustrate with a numerical example. 


2 Issue and literature review 


The DEA methodology is one of the most widely used decision aiding ap- 
proaches in the literature. Emrouznejad and Yang (2017) [5] listed more 
than 10,300 published articles. The fields of application are very diverse and 
broad. “Resource allocation is the generic problem of assigning available re- 
sources to users in the best possible way. Usually the resources are limited, 
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thus sets of activities compete for these resources establishing explicit and 
implicit dependencies, which may be subject to uncertainty. The resulting 
decision problem on the best use of the resources can be decomposed in sev- 
eral problems that have been intensively studied in the field of operations 
research in the last few decades. Examples of such problems are !: 
Scheduling problem. Scheduling is used to optimally allocate scarce resources 
to activities and includes problems like the allocation of plant and machinery 
resources, the organization of human resources, and the control of production 
processes and materials purchase. 

Knapsack problem. The main idea is to consider the capacity of the knapsack 
as the available amount of resource and the item types as activities to which 
this resource can be allocated. 

Cutting stock problem. It is a general resource allocation problem, where 
the objective is to cut pieces of stock material into pieces of specified sizes, 
minimizing the wasted material. 

Power generation scheduling problem. Power generation scheduling is re- 
quired in order to find the optimum allocation of energy such that the an- 
nual operating cost of a power system is minimized, or the obtained profit is 
maximized.” 

There is a large literature dealing with these problems and several solving 
techniques have been proposed. Several authors have been interested in the 
problem of resource allocation. Fang, [6] used a generalized DEA model for 
centralized resource allocation to uncover the sources of such total input con- 
traction. Wu et al. [14] integrated DEA and MOLP to deal with the resource 
allocation problem. Several authors have been interested in the linkage of 
the two methods, the DEA method coupled with the MOILP method (rkc et 
al. [10] and Sueyoshi [12]). We looked at questions such as: how can we take 
into account the performance indices of the different units of a production 
system to reallocate additional resources. To answer this question, we have 
made certain assumptions. First, the observed efficient production boundary 
of the system provides the best production function of the moment. Second, 
the lower expected returns to scale of a unit is an indicator to minimize the 
risk of resource waste. Indeed, when a unit has a variable returns to scale, 
the methods give a forecast of the expected yield. This random aspect led 
us to focus on the lowest expected return. In developing countries, the prob- 
lem of allocating additional resources concerns more than one. Donors and 
technical and financial partners, not only do they want to have measurable 
performance indicators, but they also want to know if an efficient unit is not 
going to lose performance if we add too many resources to it. With respect to 
inefficient units, the same questions arise. Our approach is to address these 
two types of concerns: to take into account the performance indices calculated 
by DEA and to maximize the expected returns to scale of the units taking 
into account the risks of waste. For example, if the minimum return to scale 
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is zero, then the risk of waste is high. As a precautionary measure, we have 
to wait for another evaluation period before adding additional resources. 


3 Modeling of the problem 


The Modeling of the problem requires intermediate steps to calculate, on the 
one hand, coefficients of improvement of efficiency, according to the increase 
inputs, and on the other hand, minimum values of variable returns to scale 
of the different units. 


3.1 Determination of reference units and estimation of 
the “technical” function for improving returns 


One of the major information in the process of modeling the problem is 
knowledge of reference units that are technically effective units of system. 
These units are determined using the Charnes, Cooper and Rhodes (CCR) 
model [4]. We note 


E€={j € {1,...,Ng}such as Decision Making Unit (DMU) : j is CCR efficient} , 


which contains the technical efficient DMUs. To obtain €, we use the DEA 
techniques [1,4]. Once the set € is determined, we use the Banker, Charnes 
and Cooper (BCC) model [1] with input orientation analyzing the technical 
efficient units by taking turns, individually each output factor associated with 
all input factors. This step consists of measuring the performance of technical 
efficient units relative to each output factor obtained with input factors. It 
is a question of measuring the performance of the Ng units in € with a 
single output k and the m input factors. Thus, we estimate the possible 
maximum efficiency for each output factor in relation to the observed data. 
This mathematically translates into resolution for each unit d € € and for 
each output factor k € {1,...,s} the following problem. 


max ft(k,d) = ueypa + ut 


ii CN kd kd . 
MEYkj — Yj Cig tu” <0, jg=1,...,n, 
i=l 
. . m 
: kd 
subject to the constraints (s.c) Sag = i, 

i=l 
> & ee >e, foralli,u®? eR, 


(1) 
Remark: Although the units under evaluation are all technically effective, 
therefore using the full BCC model, they are not necessarily efficient here 
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since only one output factor is considered and the unit may be less efficient for 
this single factor. In other words, the optimal value simply checks ft(k,d) < 
il 

To avoid problems related to situations of multiple optimal solutions, we 
use the following min max model (2). 


min Ad 
m 
kd kd : 
UE (Yk — Yea) — Dv; (stag — tia) — 67° < 0,5 =1,...,0, 
i=l 
m 
k,d 
s.c. ¢ 2,4; tia = 1, (2) 
i=l 
d k,d 
AY ~ 6. 2 0, jJF4, »r, 
p> se, ye Se tS Tet iag ts 


The above formulation indicates that for a given output factor k, the hyper- 
plane representing the maximum efficiency of that factor with the reference 
unit d is defined by 


m 

i; ~k,d 

[ih (Yr — Yka) — S- Be" (a = ta) = AG. 
i=1 


It can therefore be seen that, for the output factor k and for the DMU d, 
i (Ye — Yea) is provided by poor performance 


The contribution of the various input factors to this increase can therefore 


m 
5 ~k,d . is 
be measured by the expression )> ;""x;, where x; represents the increase in 


i=1 
factor input 2. In order to take into account all the reference units d € E, we 
consider the average of these values by introducing the quantity 


m 

~k,d 
y aka;, where a¥ = me. (3) 
w=1 


The coefficients ak > 0 reflect the relationship between the increase of the 
input factor 7 and the maximum average return to which it contributes for 
the improvement of the output k. 


3.2 Determination of returns to scale values 


However, it is also necessary to take into account the information on the vari- 
able returns to scale of the different units in order to allocate the additional 
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amount of available input optimally. In order to do this, we determine the 
minimum values of the returns to scale that each unit can present. These 
values are calculated for projections obtained on the efficient frontier, that 
is, if the unit is not efficient, we project it on the efficient border before cal- 
culating the minimum value of return to scale it can present. We use the 
input-oriented BCC model to determine the values of the projections and 
the minimum values of the returns to scale. We note X; and Y; the vectors 
of the inputs and outputs, respectively, corresponding to the projection of 
unit 7 on the efficient frontier. For each unit d, d=1,...,n, we calculate 
p™-4 the finite minimum value of the present returns to scale by solving the 
following problem: 


min,d d 


U = min u 


m 
yey ye oy +e? <0, jJHl,.-.,7, j#d, 
j=1 


rot I= 
“— » Gra +u?= 1, 
pe > 0; vi >0, forallr,i, u4 ER, 


According to Banker and Thrall [1], 


min,d = if 0 


p 1-— ymin,d 2 


Note y; = pind, 4=1,...,n. We note that the objective is to determine 
the smallest theoretically observable scale performance. 


4 Mathematical model 


We now describe the formulation of the model, starting with the constraints 
and then the objective functions. 


4.1 Constraints on resource availability 


We consider the following: 

—that for each input factor i, a b; quantity is available to be distributed 
between the different units 7. 

—that each DMU j expresses a request a;; units for input i (we assume a,; < 
b;, for all j. ) by introducing binary variables t;,; by t;,; = 1, if the request 
for input 7 of the DMU 9 is retained and t;,; = 0, otherwise. The constraints 
of the MOILP-DEA model are then written 
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S aejtig S0:, i=1,...,m. (5) 
g=1 


4.2 Criteria 


In the multicriteria model, we introduce s objective functions z,,k =1,...,s 
corresponding to the s output factors. z, will translate the additional return 
(expected) of the factor k due to the increase of inputs i granted to the 
different DMUs j. In order to calculate this increase it is necessary to take 
into account both the coefficients a determined in 3.1 and the coefficients 
yj; determined in 3.2. 

Also the z; function is written 


4.3 MOILP-DEA 


The multicriteria problem that we propose to allocate additional resources 
to DMUs is thus finally formulated by 


m n 


“max” zp = >> che ts5, BS D8 og 8 
i=15=1 
n 
. 7 
>) Geis =: 41s, 'M; ( ) 
S.C. j=l 


ti; € {0,1} for alli, for all j. 


It should be noted that, since no interaction is involved in the constraints 
between the variables relating to two different inputs, the problem can be 
broken down into separate problems, considering each input separately. 

For a fixed input 2, the problem can be written 


n 


“max” zhi = >> Critigs b= 1.8 
j=l 
3 <b (8) 
de aigti,g S di, 
S.C. j=l 


ti; € {0,1} for all 7. 
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These m problems are each presented as a multicriteria “Knapsack prob- 
lem” for which many methods of resolution exist. We have of course 


m 
Zk = ) Zkyi- 


i=l 


5 Illustration 


We consider the illustrative example with 3 input and 3 output factors, the 
data of which are presented in Table 1. 


Table 1: Data and Results of the Example 


DMU | Input 1 Input 2 Input 3 Output 1 Output 2 Output 3 | score CCR | p™™4 
1 8 4 6 5 4 8 1.0000 ‘| 0.6748 
2 7 8 4 3 7 4 0.9960 ‘| 0.6772 
3 11 6 4 11 11 5 1.0000 | 0.3423 
4 10 10 2 9 14 1 1.0000 | 0.0732 
5 2 7 3 4 iL 4 1.0000 | 0.9582 
6 6 10 10 8 10 13 1.0000 | 0.0000 
i 11 7 11 11 2 16 1.0000 | 0.0000 
8 5 14 8 6 13 9 0.8885 0.2899 
9 10 11 9 9 5 12 0.9101 0.8969 
10 1 8 2 9 13 2 1.0000 | 0.0388 
11 6 5 3 12 5 5 1.0000 | 0.1574 
12 5 9 1 10 8 1 1.0000 | 0.0845 
13 9 12 8 14 17 10 1.0000 | 0.0000 
14 13 13 7 15 12 9 0.8099 | 0.0000 


Using the CCR model, we get all the technically efficient DMUs 
E = {1,...,14} \ {2,8,9, 14}. 


Using the approach proposed by Banker and Thrall in [1] (see problem (4) 
above), we obtain the minimum values of the variable returns to scale p™™4 
of the 14 DMUs considered (see the last column of Table 1). To analyze each 
d € € and each output factor k, we use the problem (2) and we obtain the 
following results (see Table 2). The results in Table 3 allow the calculation 
of the coefficients a’. 


Let us consider that the requests a;; expressed in units 7 for the different 
DMUs j are those shown in Table 4. 
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Table 3: Results of the coefficients ak 


ae Trae — ae 
0.0233 | 0.0876 | 0.1560 
0.0812 | 0.0973 
0.0406 | 0.0698 | 0.1513 


wir) e] 
= 
So 
oo 
S 
a 


Table 4: Quantities requested aiz 


DMU 123 45 67 8 9 10 11 12 18 14)8; 
Inputs 7 

1 1203011021 1 0 2 #40 | 10 
2 1ii1oe3i1t10i1ti1d1 «0 2 «0 1575 
3 2001211100 2 1 1 £2 ~'58 


With this example, the MOILP-DEA model is equivalent to solving the 
3 following multicriteria “knapsack problems” . 


14 
en Sui 3 
Zia = Oy Dd) aj 71,3 
j=l 


14 
9 as, Bd, 
max” ¢ 22,1 =a] >) a1j7jt1, 
j=l 


14 
53 
Z31 =O} D) a1; 7j¢1,5 
jal 
14 


YS ar jtiy < 10, 
j=l 
#14 €{0,1} j=1,...,14 


14 
— 1 
Z1,2 = 3D) a2;7;t2,; 
j=l 
14 
“max” 4 Z2,2 = a5 D> daj7;¢2,; 
j=l 


14 
— 743 
Z3,2 = 5 )) a2;7;¢2,; 
j=l 
14 


YS agjtej <5, 
j=l 


fo € {01} GH Tyce 1A. 
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14 
— 1 
Z1,3 = 03 )) 43;7;t3,j 
j=l 
14 
“max” ¢ Z2,3 = 03 D> a3;7;t3,; 
j=l 
7 14 
Z3,3 = 03 >) a3;7;¢s,; 
j=l 
14 
a3 it3 < 8 
S.C. ae food ; 


jJ= 
tz; € {0,1} PH Aye ee DA. 


Individual maximizations of the 3 criteria give: 
k 1, t141 t1,2 tha t19 t110 ti 1, and 13 = 0 otherwise 
k 2, to1 t2,2 t2.3 to.5 t2,9 1, and t25 = 0 otherwise 
k 1, 31 t3,5 t3.8 3.11 I. and t3,5 = 0 otherwise. 


6 Conclusion 


In this article, using the relationship between DEA efficiency and MOILP, we 
have proposed a mathematical model, a multicriteria linear problem for ad- 
ditional resources allocation, and we have solved an illustrative sample. The 
formulation described here can provide an efficient framework for optimizing 
investment taking into account analyses of the performance efficiency of a 
production system. On the one hand, the determination of the units of refer- 
ence can facilitate the planning of the development programs of a given unit 
(efficient or inefficient); on the other hand, the minimum values of variable 
returns to scale make it possible to analyze the minimum need for additional 
input factors to improve returns. This different information provided by the 
use of DEA makes it possible to efficiently allocate additional resources so, 
to search and performance efficiency, and stability of efficiency of the units 
of the system. We then hope, for validation purposes, to apply these models 
to specific cases relating to public sectors such as education. In the case of 
aid and financing projects for farmers, it is important to note that this is the 
case in the case of health and a private sector such as agriculture. This would 
then allow validation of this initial model proposal and could be re-designed 
if needed. 


Acknowledgements 


Authors are grateful to there anonymous referees and editor for their con- 
structive comments. 


204 J.W. Yougbare 


References 


1. Banker R.D. and Thrall R.M. Estimation of returns to scale using data 
envelopment analysis, Eur. J. Oper. Res. 62 (1992), 74-84. 


2. Belton V. , Vickers S. Demystifying DEA- a visual interactive approach 
based on multiple criteria analysis, J. Oper. Res. Soc. 44 (1993), 883-896. 


3. Bouyssou D. Using DEA as a tool for MCDM: some remarks, J. Oper. 
Res. Soc. 50 (1999), 974-978. 


4. Charnes A., Cooper W.W. and Rhodes E. Measuring the efficiency of 
decision making units, Eur. J. Oper. Res. 2 (1978), 429-444. 


5. Emrouznejad A. and Yang G., A survey and analysis of the first 40 years 
of scholarly literature in DEA: 19782016, Socio-Econ. Plan. Sci. 61 (2018), 
4-8. 


6. Fang L., A generalized DEA model for centralized resource allocation, Eur. 
J. Oper. Res. 228(2) (2013), 405-412. 


7. Jahanshahloo G. R., Lotfi F. H., Shoja N. and Tohidi G. A method for 
finding efficient DMUs in DEA using 0-1 linear programming, Appl. Math. 
Comput. 159 (2004), 37-45. 


8. Joro T., Korhonen P. and Wallenius J. Structural comparison of data en- 
velopment analysis and multiple objective linear programming, Manage- 


ment Science. 44 (1998), no. 7, 962-970. 


9. Korhonen, K. and Syrjnen, M. Resource allocation based on efficiency 
analysis, Management Science. 50 (2004) (8),1134-1144. 


10. rkc H., nsal M.G. and Bal H. A modification of a mixed integer linear 
programming (MILP) model to avoid the computational complexity, Ann. 
Oper. Res. 235(1) (2015), 599-623. 


11. Stewart T.J. Relationships between data envelopment analysis and mul- 
tiple criteria analysis, J. Oper. Res. Soc. 47 (1996), no. 5, 654-665. 


12. Sueyoshi T. Mixed integer programming approach of extended DEA- 
discriminant analysis, Eur. J. Oper. Res. 152 (1) (2004), 45-55. 


13. Thanassoulis E. and Dyson R.G. Estimating preferred target input-output 
levels using data envelopment analysis, Eur. J. Oper. Res. 56 (1992), 80— 
97. 


14. Wu J., Zhu Q., An Q., Chu J., and Ji X. Resource allocation based on 
context-dependent data envelopment analysis and a multi-objective linear 
programming approach, Comput. Ind. Eng.101 (2016), 81-90. 


Multicriteria allocation model of additional resources ... 205 


15. Yougbar J. W. and Teghem J. Relationships between Pareto optimality in 
multi-objective 0-1 linear programming and DEA efficiency, Eur. J. Oper. 
Res. 183 (2007), 608-617. 


16. Yun Y.B., Nakayama H., Tanino T. and Arakawa M. Generation of effi- 
cient frontiers in multi-objective optimization problems by generalized data 
envelopment analysis, Eur. J. Oper. Res. 129 (2001), 586-595. 


